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Abstract. In the present paper we establish the necessary and sufficient conditions for two 
ordinary differential equations of the form y" 2 + A(x, y, y')y" + B(x, y, y') = to be equivalent 
under the action of the pseudogroup of contact transformations. These conditions are formulated 
' in terms of integrals of some one-dimensional distributions. 
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1 Introduction 

In the present paper we solve the problem of equivalence of 2nd order ODEs of the form 



o 

25: y" 2 + A(x,y,y')y" + B(x,y,y') = (1) 



under the action of the pseudogroup of contact transformations of the space M 3 (x,y,p), where 
p = y' ■ We use the results of A. Tresse and B. Kruglikov who described the algebra of differential 
invariants of 2nd order ODEs of the form 



00 

m! y" = f{x,y,y') 

O 

: under the action of the pseudogroup of point transformations of the space M. 2 (x,y). 

S. Lie showed that any two second order ODEs of the form (2) are equivalent under the 
action of contact transformations pseudogroup of the space R 3 (x,y,p). The problem of point 
equivalence of such ODEs was the point of investigation of many papers. Lie proved that if the 
ODE (2) is point equivalent to the linear equation y" = 0, it is necessary cubic with respect to 
$_i . the first derivative, i.e., is of the form 

y" = a 3 (x, y)y' 3 + a 2 (x, y)y' 2 + a x (x, y)y' + a (x, y). 

The class of such equations is closed under the action of point transformations pseudogroup. 
Moreover, he formulated the sufficient conditions of the ODE to be linearizable (see [5]). R. Liou- 
ville [6] found precise conditions for linearization. A. Tresse [9] found the complete set of relative 
differential invariants of ODEs (2). B. Kruglikov [4] formulated the results of Tresse in a modern 
language. He described the algebra of absolute differential invariants and solved the problem 
of equivalence of such ODEs under the action of point transformations pseudogroup. The class 
of ODEs, cubic with respect to the first derivatives was considered by Tresse [8], who also 
found the relative invariants of these ODEs. V. Yumaguzhin [10, 11] studied the the problem of 
point equivalence of such ODEs and the obstructions to linearization. We also note the book of 
N. Ibragimov [2], the paper of N. Ibragimov and F. Magri [3] and the paper of O. Morozov [7]. 
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2 Preliminaries 



Let 7r be a vector bundle and G be the pseudo-group of diffeomorphisms acting on ir. The action 
of G naturally prolongs to the action on the space J k ir of fc-jets of the sections of it. The function 
I G C°°(J fc 7r) is called an (absolute) scalar differential invariant of order k if it is constant along 
the orbits of prolonged action of G on J k ir. The set of all absolute differential invariants is an 
algebra A with respect to the standard arithmetic operations. 

The function F G C°°(J fc 7r) is called a relative scalar differential invariant of order k if for 
every g 6 G one has g*F = n(g) • F, for a smooth function /i : G — > C°°(J k 7r) satisfying the 
conditions fi(g • h) = h*fi(g) • fi(h), fi(e) = 1. In other words, the equation F = is invariant 
under the action of G. The function /x is called a weight. Let M = {/J,} be the space of all weights. 
Denote by W 1 the space of relative differential invariants of weight [i. Then 1Z = |J W 1 is a 

.M-graded module over the algebra of absolute differential invariants. An invariant differential 
operator is a combination of total derivatives T> x % in base coordinates (x 1 , . . . , x n ) of it 

n 

A = Y, AiV x i + A, A i} A G C°°(JV), 

i=l 

which commutes with the action of G. 

Let us formulate the classification of Tresse and Kruglikov which we will need in what follows. 

Any ODE y" = f(x,y,y') may be represented as a section of the bundle J°M. 3 (x , y , p) = 
M. 4 (x,y,p,q), where p = y', q = y" denote the corresponding coordinates in jet space. Let V x , 
V y , V p be the total derivatives in x, y and p respectively. We also denote V x = V x +pV y . Let 
lL = V l x V™V k p (q). 

The weights form a two-dimensional lattice and the space of relative invariants is a direct 
sum 1Z = rsGZ 7£ r ' s (see [4]). First relative invariants appear in order 4, they are 

/ = 4,, H = qi - Aql, + 6q° 02 + q(2qf - 3q 2 m ) - 

~ 9oo(?io - 4 9oi) + <?oo<??o - ^oo^i +Q-Q- Qoo- (3) 

Invariant differential operators are 

A p = V p + (r-s)^:TZ r ' s ^TZ r -^ s+1 , 
A x = V x + pV y + qVp + r ( 3<4 + 2 g " g ± gfo ) + 
+ a + g ° 5 ° g 4 +g M : 1W ft'+M 

y ~ Ho x + V + P Ho) y+ V %0+ < ) p+ 

r /3goo ■ goi 1 19 googoo 1 21(g 5 og + goi)g 5 o A 
+ 4g 4 + 10g& + 204,4, / 



+ s 



v 4 + 24 + 54 + 1044 ; • • 



We say that the equation y" = f(x, y, y') is generic, if /i? / for it. 
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Theorem 1. [4] The space of relative differential invariants of a generic equation y" = 
f(x,y,y') is generated by an invariant H and the invariant differential operators A x , A y , A p . 

We also formulate the equivalence theorem for the generic ODEs. 
The functions 

Hio = A X H, Hqi = A y H, K = A p H 

form the basis of invariants of order 5. There are 11 basic invariants in order 6, we need 8 out 
of them, namely 

H 20 = A 2 X H, H 11 = A x A y H, H 02 = A 2 y H, K 10 = A X K, K 01 = A y K, 

n 6 = 5&j _6.^0 j nl Q = lL {[ A p ,A x ]H-A y H), ni = AlH-^H. 

We also denote 

J, = /-l/8^3/8 G n l, 0) j 2 = jl/4 H l/4 € K 0,l 

Using J i and J 2 construct the following invariants of pure order 5 and 6: 

H 10 = H 10 /(JfJ 2 ), H 01 = H 01 /(J 2 4), K = K/(J 1 J 2 2 ), 
H 2 o = H 2 q/(J 1 J 2 ), Hu = Hn/(J 1 J 2 ), H 02 = H 02 /(J 2 J 2 ), 
K w = K 10 /(J 2 J 2 ), K 01 = 7C i/(JiJ 2 3 ), 

fi 6 = n 6 /(j 1 4 j|), = ^10/(^1 2 Ji)> ^20 = ^20/(^2 )• 

Any second order ODE £ may be considered as a section ss of a bundle J°M 3 (x, y,p). Hence, 
any differential invariant I of order k may be restricted to the equation £ via pull-back of the 
fcth prolongation of the section S£\ 

I £ := («?>)'(!). 

Let £ be a generic second order ODE such that the functions Sf , H^, K £ are local coor- 
dinates on M. 3 (x,y,p). The other differential invariants can be expressed as functions of these 
three: 

H% = ^{Hfo, H £ Q1 ,K% K% = *g-(fff , H^,K e ), fi*f = T*f(77f , H^,K £ ). 

Theorem 2. [4] Two generic second order ODEs £\ and £ 2 are equivalent under the action 
of the pseudogroup of point transformations if and only if the functions 

cb £ cb £ <b £ \Sf £ \l> £ T %£ T 5£ T 4£ 

^20' ^11) ^02' ^10' W 01> 1 i i 10' i 20 

coincide. 

The criteria for equivalence of non-generic ODEs can be found in [4]. 



3 Various formulations of the problem 

Let M be a 3-dimensional contact manifold with the contact distribution C. Let (x,y,p) denote 
the canonical coordinates in M. Then C is determined by a contact 1-form 00 = dy — pdx. 
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We will say the the ODE (1) is of hyperbolic type if D = A 2 - 4B > 0, and that it is of of 
elliptic type if D < 0. Since we are interested in the real case, we restrict ourselves to the ODEs 
of hyperbolic type. Let 

y" = Xi(x,y,y') and y" = X 2 (x, y, y') (5) 

be the roots of this ODE. Note that D is a relative contact invariant on the set of ODEs (1). 

Geometrically, each of the ODEs y" = Xi(x,y,y') determines a 1-dimensional distribution 
T% in M. 3 (x, y,p) = J°M. 2 (x,y). This distribution lies in C and is spanned by a vector field 

^ = ^+p| + A,(x, y ,p)|. (6) 

The Lie bracket [Xi,X2] does not belong to C. 

Definition. We will say that two 1-dimensional distributions T\ and T<i in a 3-dimensional 
manifold M are in general position if 

1) vector fields X, Y which define these distributions, are linear independent at any point; 

2) their Lie bracket [X, Y] is linear independent with X, Y. 

Note that the 2-dimensional distribution T\ © T% determines the contact structure on M. 

Theorem 3. The ODE (1) of hyperbolic type is equivalent to a pair of 1-dimensional distri- 
butions in general position in IR 3 (x, y,p). 
Proof. The necessity is obvious. 

Let T\ and T2 be two 1-dimensional distributions in general position in M. 3 (x,y,p) and let 
X, Y be vector fields determining these distributions. Denote T = T\ © Ti- Let us prove that 
the condition 2) of the definition does not depend on a choice of the basis in T . For any functions 
/ and g we have 

[fX,gY] = fg-[X,Y] + f -Xg-Y -g-Yf -X = fg-[X,Y] mod T. 

It follows that for another basis X' = aX + /3Y, Y' = jX + SY one has 

[X', Y'\ = (aS - £7) • [X, Y] / mod J", 

since [X, Y] ^ mod J 7 , aS - £7 ^ 0. 

Let a, b be independent integrals of X and /, g be independent integrals of Y. Since X and 
Y are linear independent, any three of these functions are independent and the fourth one can 
be expressed in terms of these three. Let, for example, g = h(a, b, /). Choose the functions a, b, 
f as local coordinates. In these coordinates the fields X and Y (up to the factors) have the form 

v 9 v u 9 u 9 

x = Tf Y = hb d- a - ha db- 

Note that h a 7^ 0, since b, f, g are independent. Similarly, hf, 7^ 0. Denote c = —h a /hb. Since 
does not belong to J 7 , one has 

h a hf - h b h a f / 0. 
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This is equivalent to the inequality 



which means that one can choose a, 
have the form (up to the factors) 

X-— Y ~ 9 1 c 9 Flex W 9 
dc da db hf dc 

where we denoted 

Flex(/i) = h aa hb ~ 2h a h b h ab + h bb h 2 a . 

The distribution T is determined by the 1-form to = db — cda and defines the contact structure 
on R 3 . 

Note that if a vector field Z lies in T and is linear independent with X, then up to a factor it 
is of the form (6). Consider the contact transformation which satisfies the following conditions: 
1) it preserves some point in M 3 , 2) its differential is sufficiently close to identity and 3) the 
direction of the field X is not an eigendirection. Such a transformation maps X and Y to the 
fields which are not proportional to X. These vector fields up to the factors will have the form 
(6). Thus, they determine two equations (5). □ 

Let us also show that the ODE (1) can be treated as two pairs of functions on M, satisfying 
some additional conditions. 

The equation y" = \\{x,y,y') locally has a pair of integrals (a, b) (which are also integrals 
of Xi) such that da A db / 0. Since X\ lies in the kernel of the forms oj, da and db, one has 
u> = a da + /3 db for some functions a, f3. Thus, ODE (1) of hyperbolic type, determines two 
pairs of functions (a, b) and (/, g) such that any three of the are independent. It follows from the 
equality da(X\) = X\a that 1-forms (X±g) da — (X\f) db (Xjfe) df — (^a) dg are proportional 
to ui. 

Definition. Let (a, b) and (/, g) be two pairs of functions on M. We will say that they are 
in general position, if: 

1) any three of them are independent; 

2) the 1-dimensional distributions determined by X and Y are in general position. 

Proposition 1. Let two pairs of functions (a,b) and (f,g) on M be in general position. Let 
a' = a'(a,b), b' = b'(a,b) be independent functions, and f = f'(f,g), g' = g'(f,g) be independent 
functions. Then the pairs (a',b') and (f',g') also are in general position. 

Proof. Follows from the fact that both pairs of functions (a, b) and (a', b') determine the 
same distributions. □ 

Theorem 4. The ODE (1) of hyperbolic type is equivalent to two pairs of functions (a,b), 
(/, g) on M in general position. 



b, c as coordinates on M 3 . In these coordinates X, Y will 
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4 The main result 



Now we return to the ODEs y" = Xi(x,y,y') and vector fields X{, i = 1,2. Choose two pairs of 
independent integrals (a, b) and (/, g) of X\ and X 2 respectively. Let g = h(a, b, /). Let 

a, b, c= - — 
h b 

be new coordinates in M 3 and let G 2 (a, b, c) = — F\ex(h)/hf. Then Xi and X 2 will have the form 

X 1 = l and X 2 = | + c| + G2 (a, 6 ,c)|, (7) 

respectively. We will say that (a, 6, c) are canonical coordinates and that the ODE 

b" = G 2 (a,b,U), 6 = 6(a), (8) 

determined by X2 is associated with the ODE (1). 

Let M3 and M3 be two 3-dimensional contact manifolds. Consider the contact transformation 
ip : M3 —7- M3 which maps the ODE (1) into another equation of the same form. Let it map the 
vector field X\ to the vector field X\ = ip*X\. We prove that <p acts of associated equations as 
a point transformation. 

Let M2 = M3/X1 be the factor of M3 by the trajectories of X\ (locally it can be done). Let 
7r : M3 —7- M2 be the canonical projection. Similarly, let tt : M3 — > M 2 = M^jX\. The modelling 
bundle for tt is the bundle 7To : J^iV — > J°N where iV is some 1-dimensional manifold. 

Let ipo : — > J l N be the restriction on tp to J l N . The diffeomorphism cp maps the 
integrals a and b of the distribution determined by X\ to some integrals a and b of the distribution 
determined by Xi, i.e., to functions of a and b. Thus, it induces the diffeomorphism if)Q : J°iV — )■ 
J°iV, which closes the diagram 

jiAT ^ ^Jijv 

JU N to ^J0 N , 

Since one can choose the functions a and b as coordinates on J°N, one can treat ^0 as a 
point transformation. Conversely, any diffeomorphism ^0 : — > J°N lifts to the contact 
diffeomorphism cpo : — > J X N. Hence, the canonical coordinate system is defined up to a 
change of (a, b) to the pair of independent functions (a(a, b), b(a, b)). 

Proposition 2. The point equivalence class of the associated equation (8) does not depend 
on the choice of the integrals (a,b), (f,g)- 

Proof. The choice of another integrals a, b of X\ instead of a, b induces the change of 
canonical coordinates, thus does not change the equivalence class of the ODE (8). 

Let F(a, b, f,g) = denote some functional dependency between these integrals. Any other 
dependency (in particular, g — h(a, b, f) = 0) has the form 

F(a,b,f,9) = *(F(a,b,f,9))=0. 
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Clearly, the ratio F a jF b does not depend on $, in particular, 

Fa = K 

F b h h 

The choice of another integrals / = 4>(f,g) g = ip(f,g) also does not change the ratio F a /Ff,, 
hence, does not change the canonical coordinate system. □ 

In the construction above we can change the order of X\ and Xi- Thus we obtain another 
ODE b" = Gi(a,b,b'), associated with (1). We say that these ODEs are dual to each other 
(see [1]). 

From above it follows the 

Theorem 5. Let £ and £ be two ODEs (1) of hyperbolic type and let ([£?], an< ^ 
[£%]) be the pairs of point equivalence classes of their associated ODEs. The ODEs £ 
and £ are contact equivalent if and only if one of the classes ([£f], [£%]) coincides with one of 
the classes ([£f], [£§]). 

Let us present some examples. 
Example 1. Consider the equation 

y" 2 - y" = o. (9) 

Here Ai = 1, A2 = 0. Let 

a = p-x, b=^x 2 -px + y, f = p, g = y-px. 

Then g = b — ^(a — f) 2 , which yields G2 = 1. On the other hand, b = g + ^(a — f) 2 , which yields 
Gi = —1. Thus, the pair of ODEs associated with (9) are 

b" = 1 and b" = -1. 

Clearly, they are point equivalent. 

Example 2. Consider the equation 

(y" - y)(y" - y ) = o. 

Here we have Ai = y, A2 = p- Let 

a = p 2 -x 2 , b = x-\n(p + y), f = y 
The pair of associated ODEs is 

y - = - y < 2 + °"'' and V 
a 

Both of them are point equivalent to the ODE b" = 0. 
Example 3. Consider the equation 

(y" + x)(y"-l) = 0, 



-p, g = pe 



a 



7 



for which Ai = — x, A2 = 1. Let 



a = ^x 2 +p, b=^x 3 -px + y, f = p-x, g = ^x 2 -px + y. 



The pair of associated ODEs is 



1 , u „ -l-46 ,2 + 2V6' 2 + 46' 4 

b = — and b 



b'-l (l + 86 /2 -4V6' 2 + 46' 4 ) 3 / 2 ' 

Both of them admit the same Lie algebra, hence they are point equivalent. 

Example 4. We give the example of two pairs of functions in general position for which 
associated ODEs are not point equivalent to each other. Let functions (a, b) and (/, g) satisfy 
the conditions 

a 2 - 2ab + f 2 + 26/ - g 2 + b 2 = 1, g > 0, b - a + / > 0. 

Then g = h(a, 6, /) = \J a 2 — lab + f 2 + 2b f + b 2 — 1, and the associated equation is 

6' 2b' 2 b' 3 
b" = r r + 



a — 6 a — 6 a — 6 

It is point equivalent to the equation b" = 0. 

On the other hand, b = h(f, g,a) = a — f + yl — 2a f + g 2 , and the second associated 
equation is not cubic in the first derivative, hence it is not point equivalent to the first one. 

In conclusion we give the formulas needed to calculate the invariants used in Tresse-Kruglikov 
d d d 

theorem. We denote — , — and — by V a , V5 and V c . In the initial coordinates (x,y,p) they 
oa ob oc 

are of the form 
where 

Z\ — CLpb x fy Q-pbyfx Q-yPxfp ~\~ Qybpfx ~i~ &xbyfp O'xbpfy, 

and the partial derivatives of / in x, y and p are understood as derivatives of a composite 
function. 

The invariant H is 

7 / h \ 2 h 

H = V 2 V 2 F - 2 • ■ V 2 c V a V b F + (-^-j • V 2 V 2 F - 4V c V a V fe F + 4 ■ ■ V C V 2 F+ 
+ 6V 2 b F + F (^2V 3 c V a F - 2 — ■ Vj?V 6 F - 3V 2 V 6 i^) - V c F^V 2 V a F - ^ • V 2 V b F- 



v a = 




~ h vfy)§- x 


+ (&p/x 


v 6 = 






+ {axfp 


v c = 


^(( a yb P 




+ (Op&a: 



4V c V 6 f) + V 3 C F (v a F - ^ • V b F^j - 3V 2 F • V b F + F 2 • V^F. 
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Invariant differential operators are 
V 5 F 

A c = V c + (r-s)- 



5V4F' 

A a = V a -^ -V b + FV C + (3r + 2s)V c F+ 
rib 



+ (2r + s) 



'F • V^F + V 4 c V a F - j- ■ V 4 c V b F^ 
VtF 



/ 6F.V c 5 F + 5V4V a F-5-^-V4V5Fx 2 
+ ( 2V C F + ^ J + (3r + 2s)-|-+ 

,V^V fe F , V^F-V C F , V^F(F • V^F + V^V 6 F) 

+ (r + 2 S )-f-V + (19r + 6a) ° „ c + (21r + 6s)- ' 



4V^F v ; lOV^F v 7 20V4F-V4F 

where X> a , X>6, X> c — are the total derivatives corresponding to the operators V , V c . 

Acknowledgement. Author wishes to express his deep gratitude to V.V. Lychagin for the 
formulation of the problem and valuable comments. 
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